Pulsar glitches: The crust is not enough 
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Pulsar glitches are traditionally viewed as a manifestation of vortex dynamics associated with a 
neutron superfiuid reservoir confined to the inner crust of the star. In this Letter we show that the 
non-dissipative entrainment coupling between the neutron superfiuid and the nuclear lattice leads to 
a less mobile crust superfiuid, effectively reducing the moment of inertia associated with the angular 
momentum reservoir. Combining the latest observational data for prolific glitching pulsars with 
theoretical results for the crust entrainment we find that the required superfiuid reservoir exceeds 
that available in the crust. This challenges our understanding of the glitch phenomenon, and we 
discuss possible resolutions to the problem. 
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Context.- Mature neutron stars tend to be extremely 
stable rotators, in some cases rivalling the best terres- 
trial atomic clocks. However, in their adolescence they 
may behave in a less ordered fashion. In particular, 
many young neutron stars exhibit (more or less) regular 
glitches, where the observed spin rate suddenly increases 
(for a recent collection of glitch data, see [1]). These 
spin-up events tend to be followed by a slow relaxation 
towards the original spin-down rate. In some cases (no- 
tably the Crab pulsar) the glitch leads to a permanent 
change in the spin-down rate, but in most cases the sys- 
tem resumes spinning down as before the event. 

The archetypal glitching neutron star is the Vela pul- 
sar, which has (since the first observed event in 1969) 
exhibited a reqular sequence of similar size glitches. The 
consensus view is that these events are a manifestation 
of the presence of a superfiuid component in the star's 
interior [2]. This idea was first put forward by Ander- 
son & Itoh [3] who envisaged a glitch as a tug-of-war 
between the tendency of the neutron superfiuid to match 
the spindown rate of the rest of the star by expelling 
vortices and the impediment experienced by the moving 
vortices due to "pinning" to crust nuclei. Strong vortex 
pinning prevents the neutron superfiuid from spinning 
down, creating a spin lag with respect to the rest of the 
star (which is spun down electromagnetically) . This sit- 
uation cannot persist forever. The increasing spin lag 
leads to a build up in the Magnus force exerted on the 
vortices. Above some threshold pinning can no longer be 
sustained, the vortices break free and the excess angular 
momentum is transferred to the crust. This leads to the 
observed spin-up. 

Several decades have passed since the two-component 
model was first suggested, yet there has been surpris- 
ingly little progress on making the model quantitative. 
Most thinking has gone into the microphysics, especially 
concerning the interaction between neutron vortices and 



crust nuclei [4], a key ingredient in the scenario. Mean- 
while, the detailed mechanism responsible for triggering 
glitches in the first place remains unknown [5]. Moreover, 
we do not actually know the location of the superfiuid 
reservoir associated with these events. It is generally as- 
sumed that the vortices are pinned in the crust and hence 
the angular momentum available is that associated with 
the crust superfiuid. This notion is supported by data for 
frequent glitchers (Vela being the prime example). The 
analysis by Link et al [6] suggests that glitches represent a 
self-regulated process that involves a superfiuid reservoir 
with moment of inertia I n /I ~ 1% where I and / n are, 
respectively, the moments of inertia of the entire star and 
the superfiuid component. The similarity of the inferred 
I n to the theoretically estimated moment of inertia of the 
crust (which is dominated by the free neutrons in the in- 
ner crust) for realistic equations of state [7] supports the 
notion that glitches involve only the crust region. 

We argue that this logic breaks down when one ac- 
counts for the non-dissipative entrainment coupling be- 
tween the neutron superfiuid and the crust lattice, an 
effect which can be expressed in terms of an effective neu- 
tron mass, m*. Recent work indicates that this effective 
mass may be significantly larger than the bare neutron 
mass, m n [8, 9]. This implies a decreased mobility of the 
superfiuid with respect to the lattice and the need for a 
larger angular momentum reservoir for glitches. Combin- 
ing the latest data for glitching systems [1] with a general 
relativistic multifluid model that includes entrainment, 
we find that the required superfiuid moment of inertia 
is above the capacity of the crust superfiuid. This sug- 
gests that some fraction of the core superfiuid must par- 
ticipate in glitches, which raises a number of interesting 
questions. 

Phenomenology and observations.- The discussion of 
vortex mediated glitches is usually based on a "body" 
averaged model with two components. The first repre- 
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sents the charged component (including the elastic crust) 
which is spun down electromagnetically. Labelling this 
component by an index p, we have 

I p tl p = -afip - Mpin - A/" MF (1) 

where the first term on the right-hand side represents the 
standard torque due to a magnetic dipole (the coefficient 
a depends on the moment of inertia, the magnetic field 
strength and its orientation; we assume that these param- 
eters do not evolve with time). We also have a superfluid 
component, with index n, which evolves according to 

I n (l a =A/"pi n +A/"MF (2) 

On the right-hand sides of these equations we have added 
terms representing torques associated with vortex pin- 
ning (.A/pi n ) and dissipative mutual friction (A/mf) asso- 
ciated with scattering off of the vortices in the superfluid. 
We will not need explicit forms for these in the following. 

Glitches can be understood as a two-stage process. In 
the first phase the superfluid vortices are pinned. This 
means that A/" p i n is exactly such that (l n — 0. That is, 
the pinning force counteracts the friction which tries to 
bring the fluids into co-rotation. The upshot is that the 
crust evolves according to 

1 1 2cj 
7 p f2 p = -a^ _^ (<_ <0 ) ( 3 ) 

Sip S2q i p 

Assuming that a system starts out at co-rotation (with 
spin fio at time to), we can estimate how the spin-lag, 
Ail = il n — J7 p , between the two components evolves with 
time. As long as the spin-lag is small we have Asl/il p ss 
t g iitch/2r c where t g utch is the interglitch time and r c = 
— n p /2(l p is the characteristic age of the pulsar. 

At some point, this lag reaches a critical level where 
the vortices unpin. The two components then relax to 
co-rotation on the mutual friction timescale (which may 
be as fast as a few hundred rotations of the system [10]). 
This transfers angular momentum from the superfluid 
reservoir to the crust, leading to the observed glitch. As- 
suming that angular momentum is conserved in the pro- 
cess (such that the entire spin-lag Ail drives the observed 
glitch jump Ailp) we have 

JpASlp = I n ASl — ► ^p w ^W^ ( 4 ) 
p p lip I 2t c K J 

where I = I n + I p is the total moment of inertia (we have 
assumed a small superfluid reservoir, i.e. I w I p ). 

Let us compare this model to observations. To do this, 
we assume that we see a number of glitches in a given sys- 
tem during an observation campaign lasting t h s - Then 
we can work out the accumulated change in the observed 
spin due to glitches, and relate the result to the simple 
two-component model. From (4) we then have 

I n / 1 w 2t c A where A=—[J2 A ^p/^P J (5) 

tabs \ i J 



PSR 


t c (kyr) 


A (xl0 _9 /d) 


In/I (%) 


J0537-6910 


4.93 


2.40 


0.9 


B0833-45 (Vela) 


11.3 


1.91 


1.6 


J0631+1036 


43.6 


0.48 


1.5 


B1338-62 


12.1 


1.31 


1.2 


B1737-30 


20.6 


0.79 


1.2 


B1757-24 


15.5 


1.35 


1.5 


B1758-23 


58.4 


0.24 


1.0 


B1800-21 


15.8 


1.57 


1.8 


B1823-13 


21.5 


0.78 


1.2 


B1930+22 


38.8 


0.95 


2.7 


B2229+6114 


10.5 


0.63 


0.5 



TABLE I: Inferred superfluid moment of inertia fraction for 
glitching pulsar which have exhibited at least two (large) 
events of similar magnitude. The data is taken from [1] (up- 
dated to included a few more recent events [11]), c.f., Figures 1 
and 2. We give each pulsars name, the characteristic age, r c , 
the averaged rate of spin-reversal due to glitches, A, and the 
moment of inertia ratio I n /I obtained from (5). 

For systems that have exhibited at least two glitches of 
similar magnitude [1] we can estimate the average rever- 
sal in spindown due to (large) glitches per day of obser- 
vation, A. This leads to the inferred moment of inertia 
fractions listed in Table I. For some systems, like the 
Vela pulsar and the X-ray pulsar J0537-6910, the esti- 
mate should be quite reliable given the number of glitches 
exhibited and their regularity. In other cases, the data is 
less impressive, as is clear from Figure 2. Nevertheless, 
the message seems clear: Glitches require the superfluid 
component to be associated with at least 1-1.5% of the 
star's moment of inertia. This agrees with the conclu- 
sions of [6]. In addition, the data seems consistent with 
the idea of an angular momentum reservoir that is com- 
pletely exhausted in each event. If this is not the case 
then it is difficult to explain why the recurring glitches 
have such similar magnitude. 
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FIG. 1: The accumulated ^ AQ' p /Q p (xl0~ 9 ) as a function 
of Modified Julian date for the X-ray pulsar J0537-6910 and 
the Vela pulsar (B1833-45). The fits that lead to the slopes 
(A) listed in Table I are shown as straight lines. 

The role of entrainment- Let us now ask what the 
influence of a "heavy" superfluid may be. That is, let us 
account for the entrainment coupling. At the level of the 
averaged two-component model, the entrainment can be 
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FIG. 2: Same as Figure 1, for pulsars with a smaller number 
of large glitch events. 



expressed in terms of a coefficient e n . The two equations 
of motion then take the form [12] 

(ip - e„/ n ) tip + £nlnti n = -afip - Wpin ™ M M F (6) 



and 



(1 - e n ) I n £L n + e n I n £l p = Af pin + A/mf 



(7) 



Combining these (noting that the right-hand side of (7) 
vanishes for perfect pinning), we see that the crust now 
spins down according to 



m p = - a n 3 p 



where I = J n — 
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(8) 



Expressing the entrainment in terms of the (average) ef- 
fective neutron mass, we have 



1 



I 



(ml)' 



(9) 



The interpretation of this result is quite simple. The en- 
trainment encodes the mobility of the superfluid neutrons 
relative to the other component. If the effective mass is 
large, then the two components are effectively locked. 
Hence, the system spins down as one body (/ — > I) in 
the limit where (m*) 3> m n . Basically, the entrainment 
lowers the "effective" moment of inertia associated with 
the superfluid. 

In terms of the observed spin-down the entrainment 
only has the effect of altering the inferred magnetic field. 
We can still introduce the characteristic age (obtained 
from observables) to get the accumulated spin-down of 
the crust. The main difference now is that O n also 
changes (even when vortices are pinned). From (7) we 
have 



o n = - 



i 



-o p = ( 1 - 



<m*>,' " p 



(10) 



This has repercussions for the estimated glitch jumps be- 
cause the spin-lag between the two components takes a 



longer time to develop if the effective neutron mass is 
large. Working out the accumulated spin-lag and assum- 
ing angular momentum conservation during the glitch, 
we have (again assuming I ss I p ) 
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(11) 



The observations then provide us with the constraint 

(12) 



I m„ 



In other words, if the (average) effective neutron mass is 
large, then the constraint inferred from glitch observa- 
tions will be more severe than previously assumed (e.g., 
in [6]). This argument may not be new [13], but the effect 
has not previously been quantified. 

Moments of inertia.- In order to quantify the con- 
straint set by the glitch data, we need a relativistic model 
for the involved moments of inertia [14]. The need for 
such model is emphasised by the recent results of Chamel 
[8, 9], which suggest that the effective mass for the super- 
fluid neutrons that permeate the inner crust may, indeed, 
be quite large. The phenomenological (body averaged) 
entrainment model from the previous section illustrates 
how a large effective mass affects the analysis, but we 
need to connect this argument with a detailed neutron 
star model, incorporating realistic crust and core physics, 
as well as reasonable superfluid parameters. We build 
on the relativistic superfluid formalism developed in [14], 
and use the same assumptions regarding the supranu- 
clear equation of state and the singlet pairing gap for the 
neutrons as in [16]. The effective neutron mass in the 
crust is estimated using the phenomenological fit to the 
entrainment from [15]. We consider two models, based 
on [8] and [9], respectively. The main difference is that 
the effective mass peaks at lower densities in the latter 
case. This serves to weaken the effect we are discussing 
slightly. 

Defining the moment of inertia through J = 7f2, we 
get the total moment of inertia from [7] 



2/ , 



where R is the radius of the star, and 
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/o = T 



r i e (X-u)/2 ( p + p ) dr 



(13) 



(14) 



where p and p are the pressure and the energy density, 
respectively, and A and v determine the spacetime metric. 
In the case of the superfluid crust neutrons we ignore 
the effect of the rotational frame-dragging (this should 
be, at most, a 20% correction). Making contact with the 
phenomenological model discussed previously, we express 
the superfluid moment of inertia as (see [14]) 



J n I n Q n -\- £ n I n (Op O n ) 



(15) 
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Here 

*n«^ frW-W^dr (16) 

where n n is the number density of the free neutrons and 
/i n is the corresponding chemical potential, and 

£ * = T^T T r 4 e^/ 2 n n (m n - m n ) (17) 

In these integrals we only account for the crust super- 
fluid (R c represents the crust-core interface), but it is 
obviously straightforward to introduce a core component 
as well (as long as we keep in mind that the entrainment 
is rather different in the core). In the inner crust we can 
safely assume that the free neutrons are non-relativistic, 
which means that /i n w m n . It is then straightforward to 
work out the effect of the entrainment (in terms of e n , or 
the averaged ratio (m*)/m n via (9)) and check to what 
extent a given neutron star model satisfies the constraints 
set by the observations. Let us focus on the results ob- 
tained for the equations of state for core and crust used in 
[16]. For this model we find that I n / 1 ranges from around 
6% for a 1.2M Q star to 4% for a 1.4 M Q star and less 
than 2% for a 1.8M© star. This means that, as long as we 
do not worry about the entrainment, these models easily 
satisfy the constraints set by the observed systems; The 
angular momentum reservoir exceeds the requirements. 
Testing a few different core equations of state we find 
that this conclusion holds in general. It is worth noting 
that the I D /I ratio is significantly smaller for the models 
considered in [6, 7]. Turning to the entrainment effect, we 
find that the two models (representing the results from 
[8] and [9]) lead to (m*)/m n in the range 4 — 6 (more pre- 
cise estimates are not justified given that we are using a 
phenomcnological fit to the small number of data points 
given in [8, 9]). At this point the standard glitch logic is 
in trouble. Unless the observed systems are all low mass 
neutron stars, the crust is not enough; The associated 
superfluid does not have sufficient moment of inertia to 
explain the observations. 

Discussion - Our analysis casts doubt on the standard 
model of pulsar glitches being driven by a superfluid 
reservoir confined to the inner crust of the star. This 
seems inevitable provided the effective mass associated 
with entrainment is indeed as large as suggested by re- 
cent work [8, 9]. Basically, the entrainment lowers the 
effective superfluid to total moment of inertia by a fac- 
tor ^4 — 6 below that assumed in previous work [6]. 
This brings the theory into conflict with the data for 
pulsars with regular large glitches, as shown in Table I. 
In order to be "consistent" we need stellar models with 
In/ 1 6 — 10%, but this is only borderline achievable for 
proposed realistic equations of state provided that these 
neutron stars all have low mass. 

There are (at least) three possible interpretations of 
this analysis. The first possibility is, perhaps, the least 



attractive. The neutron star core is expected to con- 
tain superfluid neutrons in abundance. In the case of the 
singlet pairing gap we have considered (see [16]) the to- 
tal moment of inertia fraction would be at least a factor 
of two (up to an order of magnitude for massive stars) 
larger than required (and one should probably add the 
triplet pairing region to this) . This is where a new prob- 
lem arises. The remarkable regularity of the glitches in 
systems like Vela and PSR J0537-6910 suggests a, more 
or less, completely recycled reservoir of pinned superfluid 
- the very reason why the previous moment of inertia 
constraints [6] were taken as evidence in favour of crust- 
only based glitches. However, if the core superfluid takes 
part, then why are the glitches not larger (at least occa- 
sionally)? Could it be that the unknown glitch trigger 
mechanism is associated with some characteristic length- 
scale that leads to only a fraction of the total reservoir 
being replenished in each event? To make progress on 
these questions, we need more detailed work on plausi- 
ble glitch triggers (possibly associated with the onset of 
superfluid turbulence [5, 17]). 

The second possible explanation is that the core su- 
perfluid is involved in the glitch, but that the combined 
superfluid reservoir is just large enough to explain the 
observations. If this "fine-tuning" resolves the problem, 
then a more detailed calculation would be able to con- 
strain the singlet pairing gap for neutrons. This would be 
a very interesting complement to the recent constraints 
on the core superfluids obtained from the cooling rate of 
the young neutron star in the Cassiopeia A supernova 
remnant [18, 19]. 

The final possibility is, conceptually, the simplest. The 
discrepancy between the model and the observations may 
be due to a lack of "precision" in the theory. This is ob- 
viously testable. However, a more precise analysis of the 
problem requires a consistent crust-core equation of state 
along with the effective neutron mass for the superfluid 
component. So far the only results in this direction are 
those in [8, 9]. Work aimed at confirming (or refuting) 
those results should be strongly encouraged. 

There is, however, a problem also with this explana- 
tion. Unless the superfluid is confined to the crust, one 
would have to explain why the crust component decou- 
ples from the core during the glitch event. This would be 
particularly vexing if the singlet pairing gap is such that 
the neutron superfluid reaches far into the core (which 
is the case in the model we have considered), leading 
to interesting questions concerning the physics associ- 
ated with the transition from crust to core. A central 
issue concerns the nature of superfluid vortices extend- 
ing across this interface. The standard picture is that 
the vortices are magnetized in the core (due to entrain- 
ment and the presence of superconducting protons) [10] 
but not in the crust. This hints at a more complicated 
behaviour at the interface than is usually assumed. Per- 
haps there is an argument for why the crust component 
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decouples. This problem needs further consideration. 

Our analysis points to a number of interesting direc- 
tions for future work. There is an obvious need for more 
precise theoretical limits on /„//, both from theory and 
observation. On the theory side, the effective neutron 
mass in the crust should be a key ingredient in any model. 
From observations, continuous monitoring of glitching 
pulsars will improve the constraints on the supcrfluid 
reservoir. We also need to develop a better understanding 
of the superfluid/ vortex physics at the crust-core inter- 
face. Progress along these lines is required if we want to 
take our understanding of glitches to the next level, and 
use precision observations from the next generation of 
observational facilities to probe supcrfluid neutron star 
physics. 
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